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INTRODUCTION. 

Jl. This report was prepared by Dr. H. Bat email for publication by the National Advisory 
Committee for Aeronautics. The theory of small oscillations about a state of steady motion 
which was developed many years ago by E. J. Routh 1 has been applied with marked success 
in aerodynamics, the desired simplicity of the equations being secured by the introduction of 
the resistance coefficients by G. H. Bryan. 3 This simplification of the equations is based on 
the assumption that in a slight departure from a state of steady motion the increments in the 
component aerodynamical forces and couples can be expressed in terms of the increments of 
the component velocities of translation and rotation alone without any additional terms depend- 
ing, for instance, on the increments of the accelerations. This assumption seems to give a 
good approximation to the truth in the case of an airplane, but in the case of a balloon the 
additional terms are required. When a flying machine is treated as a rigid body the general 
type of steady motion is one in which the center of gravity describes a helix and the algebraic 
equation which determines the temporal characteristics of the oscillations is of the eighth 
degree, but this equation can be simplified in certain cases. In the case of an airplane having 
a plane of symmetry, the oscillations about a state of steady rectilinear flight can be regarded 
as built up from longitudinal and lateral oscillations which are practically independent of one 
another. When certain resistance coefficients are assumed to be zero each set of oscillations is 
associated with an algebraic equation of the fourth degree. 3 

A notable simplification also occurs in the ease of a body like a parachute which has an axis 
of symmetry, when the steady motion is rectilinear and in the direction of the axis of symmetry. 

In a recent report on the dynamical analysis, 3 Messrs. Klemin, Warner and Denkinger 
have studied the effect on the period and rate of subsidence of the pitching oscillations of an air- 
plane of a change in one of the resistance derivatives when all the others are kept constant. It 
occurred to the author that it might be worth while to continue this work by considering simul- 
taneous variations in two or three of the derivatives and to pay special attention to cases in which 
a slight change in some of the derivatives has (1) no effect, and (2) a marked effect in the charac- 
teristics of the long oscillation. 

The results which have been obtained are exhibited by means of diagrams in which two of 
the resistance coefficients are used as coordinates, and curves are drawn along which the modulus 
of decay of a long oscillation has a constant value, numbers being given to indicate the value 
of the period at various points of each curve. At Mr. Warner's suggestion numbers have been 
given in some cases to indicate the ratio of the period to the time of damping to half the initial 
amplitude, as this quantity is adopted as a measure of stability in the report to which we have 
just referred. 

i The Stabuity of a Given State of Motion (A darns Prize Essay). London, IST7. Advanced Rigid Dynamics. Chaps, m and VI. 
J Stability in Aviation. Macmiltan <Si Co. (1911). 

« Third Annual Report o£ the National Advisory Committee for Aeronautics. Tfo. 17, p. 330. (1917.) This report mil be referred to subse- 
quently as 17. 
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The method has also been adapted to the lateral oscillations of a symmetrical airplane and 
to the oscillations of a parachute. The graphical method used here is inferior in many respects 
to the beautiful one devised by L. Bairstow and J. L. Nayler, British Advisory Committee's 
Report No. 116, 1915, but the work was completed before this report came to the writer's 
notice. 

§2. Pitching oscillations. — When the pitching oscillations are regarded as independent of 
those in roll and yaw, the biquadratic equation which determines the temporal characteristics 
of the oscillations may be written in the form 

X G4\+y)(X+z)(A+w)+aj\(X+z)+X 9 y sin d^J + ^\(A\+y) 

(1) +(£c+fw)(XS+cos 8 0 ) + (tf-zx) |? (XS-sin 9 0 ) +Xf(cos 0 o -r?)=o 
where x = - V M a , y = — M q , z = — -X u , w = -Z m 

(2) V=U+Z q + X q 

and the notation is the same as in the reports of Hunsaker, Klemin, Denkinger, and Warner. 

Let X 2 — 2aX+7 be one factor of the expression on the left hand side of equation (1). 
Replacing X 2 by 2a\ — y we can reduce the above expression to one which is linear in X and 
equate to zero the terms with and without X. This gives us the two equations 

(3) y[y+A{w+z) +4Ace]= (y +2Aa)(w+2a)(z+2a) +fr(y +2Aa) 

(4) Ay 2 - y[(y +2Aa) (w +z + 2a) +A(wz + +x} + cos 8 0 (& +$w) --^ sin B 0 (tf-zx) = o, 



which generally determine % and y uniquely 1 when a and y are given, that is, when the period 

p = , ^ X ■ : and the co-efficient of subsiden 
y-y — cr 

quantity it is convenient to use the time 



p = , ^ X ■ : and the co-efficient of subsidence = a of the oscillation are given. Instead of the latter 
y-y — a 3 



Iog e 2 .09 
— a —a 



which represents the time which it takes for the amplitude of a simple oscillation to fall to 
half value. 

With the aid of equations (3) and (4) the curves t = constant (a = constant) have been 
drawn in the (x, y) plane for various values of z, w, £, ij, f, 5, and 6 0 , the value ..4 = 100 being 
adopted in each case. We can use the same diagrams for any other value of A by simply alter- 
ing the scales for x, y, and It should be noticed in fact that equation (1) is still satisfied if 
we replace A, x, y, f by kA, kx, Ky, /cf, respectively, keeping the other quantities the same. 

In diagrams I-V there are two sets of curves corresponding respectively to the values n = 1 
and 17 = 2. Each set of curves is made up of three pairs, the two curves of each pair correspond 

* They may fail to do this for certain particular values of a and y when the two linear equations in z and y are the same. It should be noticed 
however, that when h=<s, S=o, i=o, the equations give 

(2a + MI+Z)[i» s +!M(!!-l)+£(i9-l) J ! = ° 

When 2a,+w+z=o wehavsy=£i!+tf)Z and the equation of the line along which both a and7areconstantisi=o. This is one of the boundaries o( 
the region of stability in the (i, y) plane. The second factor wHwztij— I) + {(>?— II s is generally positive with the values of the resistance coeiTicient 
usually found for an airplane. If it could vanish there would be an infinite number of straight lines along which a and y are constant and connected 
by the relation 

9—1 
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to the same value of t, but the upper curve corresponds to the value f = l and the lower curve 
to the value f = o. 

The values of z, w, £, in diagram I are roughly those found for the Curtiss J. N. 2 in Report 
17. A few slight changes have been made to facilitate the calculations, but these do not affect 
the general conclusions. The values of i? in diagram I are greater than the number derived 
from the value of X w given in Report 17. This number is about 0.6 and it is remarked on page 
332 that X a decreases as the angle of incidence increases and the speed decreases. Since X a 
may actually become negative as the critical speed is approached, curves have also been drawn 
for the values tj = o and 77= —1. 

In studying these diagrams it should be noticed in the first place that the lines t= « and 
x = 0 limit the region of stability. The first of these lines is curved. When 77 is decreased this 
line rises and the region of stability becomes more restricted. On this account alone we can 
expect a decrease in 17 to be unfavorable when the other resistance derivatives remain, const ant. 
This agrees with the result found in Report 17 but it is worth while to consider the matter more 
fully. 

If we take any value of 77 and begin to draw the curves t = constant, starting from t = <*> and 
gradually decreasing t, we find that a certain minimum value t 0 is reached below which the curve 
no longer lies in the part of the region of stability shown on the diagram. It is clear from the 
diagrams that f 0 increases when 77 decreases. This increase in the time of damping is partly 
offset by an increase in period; but, if the other resistance derivatives remain constant, the ratio 
of the period to the time of damping apparently decreases with tj. This is seen more easily by 
looking at diagrams Ic and lb. 

If some of the other resistance coefficients alter at the same time as 17 the unfavorable 
effect of a decrease in the value of n may be partially or completely offset. Thus if, when we 
decrease 17 from 2 to 1, we increase y so as to keep t constant, we increase the period p and so 
improve the stability as far as the long oscillation is concerned. With small values of ij a con- 
siderable increase in y may be needed to keep t constant when 17 is decreased, but a much smaller 
increase in y may be sufficient to offset the unfavorable effect of the decrease of 77. The effect 
of decreasing 77 may also be offset by decreasing x ; for, if we decrease a; so as to keep t constant, 
the period p is seen to be increased. This may be seen very clearly in diagram I if we start with 
the point f = 13.8, p =20.1, 77 = 2, and pass first to the point with the same coordinates in the part 
of the diagram corresponding to 77 = 1 and then proceed along a line y = constant until the curve 
^ = 13.8 is reached. A comparison of the different diagrams indicates that the general form 
and arrangement of the curves is roughly the same for the different sets of values of 2, up, and £. 
An exception occurs in the case of the curves for 77 = —1 in diagram la. It will be seen that in 
this case the curves f = 0 cross the curves f =1 while in the other diagrams a curve f = 0 remains 
below the corresponding curve f = 1 . 

It is generally assumed that f = 0, since the moment about the center of gravity of the 
airplane due to air forces is zero in horizontal flight and therefore will not be affected by a 
change in speed. It is easy, however, to imagine some arrangement which will make M u and 
therefore f different from zero. The inclination of a flap held by a spring and exposed to the 
air will vary with the speed of the airplane, consequently the force on it will not be propor- 
tional to the square of the speed, and so f will be either positive or negative. 

It will be seen from the diagrams that a positive value of f is generally unfavorable to 
stability. When 7? is negative, a machine with a moderately large value of x may be an exception 
to this rule; but, when X w is positive, the stability can apparently be improved by making f 
negative. 

If this is done with the aid of a flap held by a spring, 5 changes may also be produced in 
the other resistance derivatives, particularly X a and Z u . We must therefore also study the effect 
on stability of changes in £ and z. Whis may be done with the aid of diagrams I-V. 

5 The stability of an aeroplane which has springs in the control connections has been discussed hr a more rigorous ma nner h j L. Bairstow and 
R. Jones. British Adrfsory Committee for Aeronautics. (R— M No. 210.) 

54SS9— 21 23 
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It appears that an increase in % from 2 to 3 has practically no effect on the position of the 
curves f = constant, but it does lower the period; hence an increase in £ seems to be unfavorable. 

A comparison of diagrams IV and V indicates that an increase of z from 0.1 to 0.2 lowers 
the curves t — constant and is on the whole favorable to stability. 

The general conclusion then is that, if some kind of a flap held by a spring is to be used 
to improve the stability, as far as the long oscillation is concerned, it should be chosen so as 
to decrease Z u , increase — X u , and make M u positive. 

The diagrams may also be used for other purposes; in particular, they may be used to 
confirm some of the conclusions in Report 17. It will be noticed that as a point moves toward 
the origin along a curve t = constant, the period p increases, and so the stability is improved. 8 
It may then be advantageous to decrease both VM„ and — M q , but — M q must not be decreased 
too rapidly. There is a limit to the ratio of — dy to — dx if we wish the change to be favorable 
to stability, and this limiting value may be shown clearly on the diagram by drawing the curves 

for which j is constant, in accordance with a suggestion made by Mr. Warner. These curves 

are shown in diagrams lb and Ic. It will be noticed that the curves ^ constant are steeper 

than the curves t = constant, but are not quite as steep as the curves p = constant. 

A comparison of diagrams I and II indicates that an increase of w from 2.9 to 3.9 increases 
the period p and has more effect on the time of damping t when i? is negative than when ij is 
positive. When ^ = 2, the time t apparently increases with w; but when ij = l, t increases for 
small values of t, e. g., those less than 13. S, and decreases for large values of t. When t? = 0 and 
i? = — 1, t decreases as w increases, the effect being quite marked when rj = — 1. 

In diagram VI the effect is shown of making S different from zero. It appears that by 
making B positive we decrease the time of damping and produce very little change in the 
period. Other things being equal, an increase in 5 seems to be favorable to stability. It should 
be noticed, however, that X q occurs in the expression for V, consequently unless an increase 
in X q is balanced by a decrease in Z q the numerical value of V will bo lowered and the values 
of the quantities depending on V altered. It appears then that the effect of an increase in 
X q may be offset by the change in V. It should be observed, however, that, when X q is 
about 3, the percentage change in the value of V is not large, and so the effect of the increase 
in 5 should predominate over the effect of the change in V. In diagram VI cr is written in 
place of 5. 

In diagram VII the effect is shown of a change in 6 0 , other things remaining the same. 
It appears that the stability is improved by making 0 o positive, for the curves t = constant are 
lowered. By making d 0 negative, the time of damping is increased. In particular, if the 
speed were not increased, the time of damping would be greater in a vertical dive than in 
horizontal flight. This may be seen when the stability of the parachute and helicopter are 
studied. (See Report No. 80.) 



s It should be noticed that the nearer the representative point is to the origin of coordinates! ho more marked is the effect on p and f of 
small changes in i and y. 
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Lateral oscillations. — When the lateral oscillations are analytically independent of the 
pitching osculations the biquadratic equation -which determines the temporal characteristics 
of the oscillations is 

X(X + st) (X+w) (\ + v) + tWK+x)-£y\(\ + w)+2Ck + T 1 )-&\-yt=o (1) 

where 7 

N^-faKf, N r --nZx?, N.—^KS, , (2) 
Y v = o, Y r = o, Y v =-w. 

Writing X 2 — 2aX + #=o and reducing the above equation to a linear one in X we find, on 
equating to zero the coefficients of X and 1, that 

d(w + v + x + 4a) = (o: + 2a) (w + 2a) (ij + 2a) +f(« + 2a) + 2a) + z(l -£), 

8 l -8[(x + 2a) (w + r] + 2cx)+Wr 1 +r-ty] + ZV-yt=0. (3) 

With the aid of these equations the curves t= constant may be drawn in the Cx, y) plane 
for various sets of values of z, f, tj, and 10. It will be seen from the diagrams that there is a 
marked difference between the curves for which £ is positive and those for which £ is zero or 
negative. 

It has been pointed out recently by Prof. E. B. Wilson 8 that in normal flight at a fairly 
low angle of incidence the quantity N p is negative, and that consequently there is no need 
for a fine adjustment of the values of N v , L v , N v , and L v in order to secure both spiral stability 
and stability in the Dutch roll. At larger angles of incidence N p may be zero and even positive, 
and then there is more need for a fine adjustment. On this account most of the diagrams have 
been drawn for zero and positive values of £, one diagram being deemed sufficient to indicate 
the general arrangement of the curves when £ is negative. It will be noticed that when £ is 
negative the period p is rather short but the time of damping t can also bo made short and so 
the short period is not a great 'disadvantage. By making £ positive the period can be more 
than doubled, but this is offset by the lengthening of the time of damping. 

It should be noticed that in the diagrams the region of stability is bounded by the curve 

t = <x> and the line ?/ = ~p If the representative point lies to the left of the curve t = <x> , the air- 

plane is unstable in the Dutch roll; while if the point lies above the line V = ~^> the airplane is 
spirally unstable. 9 The shape of the curves t = constant is very much the same for the different 
values of jj, f and z used in the diagrams, but the position of the line V = -j 'varies considerably. 

An examination of diagrams VIII-XII indicates that it should be possible to construct an 
airplane which, when flying at a large angle of incidence, is spirally stable and has a period for 
the Dutch roll of about 15 seconds with a fairly short time of damping. 

It has been remarked that in a moderately stable airplane the period is not sensitive, com- 
paratively large changes in the airplane having small effect on the period and damping. 10 As 
an illustration of this phenomenon it is of some interest to consider the lines along which both 
p and t are constant. Such a line is obtained by choosing a and 9 so that the two equations (3) 
are the same. Writing down the conditions 

(w + 2 a)(r} + 2a) + £-6 = - £(w + 2a) _ g(l - £) -6 (w +77 + 2a) 

. . - . — — ^ — ■ - . — . — i.— — _ „ . — , : : ^ is - ^ 

'< Cf. J. C. Hunsaker, Dynamical Stability of Aeroplanes, Smithsonian Miscellaneous Collections (Washington), Vol.62, No. 6, June (191(1), pp. 
55-57. 

8 Fourth Annual Eeport of the National Advisory Committee of Aeronautics. No. 26. Washington (1919). 
8 For these terms see Hunsaker (loc. cat.). 

10 Cf. W. S. Farren, "Full Scale Aeroplane Experiments," The Aeronautical Journal, February, 1919, p. 56. 
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we hare two equations to deter min e a and d when z. «?, £, 57 and f are given. If on the other 
hand £, f, w, a, and 0 be regarded as given, the resulting equations will determine n and 2 
uniquely. The equation for 17 is 

- ( * + 2a) (w + 2a)f = {fl> + fl[£a? (w + 2a) - If - f 1 + f*. 

If to -f 2a is small, a comparatively small percentage change in 0 means a comparatively large 
percentage change in the value of 17. For example if to = 0.1, f =0.6, £ = 2, w + 2a= — .05, 

6 = .27, 17 = .72, 2 = 1.782 
6 = .26, ■>? = .97, z=1.794 
fl = .25, 17 = 1.233, 2 = 1.777 

In diagram VIII with £ = 1 , the equation of the straight line along which p and t are constant 
is DS84r + .0241y-.3843=o; and we have u?+2a= - .0241, 0=.4953, i= 11.12, p = 8.96. 

It appears from the diagrams that the greater the curvature of the curve f = constant, the 
more rapid is the variation of p as the representative point moves along it. If the curve is 
concave to the axis of y, p increases as the point moves up the curve, while, if the curve is 
convex to the axis of y, p increases as the point moves down the curve. 

If we wish an airplane to have resistance derivatives for which the period is not sensitive, 
a good plan is to bring the representative point as close as possible to the straight line along 
which p and f are constant. When f is negative this line does not seem to appear on the 
diagram. 

A study of diagrams VIII-XII reveals the following properties : 

(1 ) The greater the value of | the greater seems to be the effect on the damping of a change 

in w. 

(2) An increase in w decreases the time of damping but does not greatly alter the period. 

(3) An increase in 17 decreases the time of damping and increases the period when £ = 0/ 
but, when £=1, 2, or 3, the effect seems to be reverse. 

(4) An increase in 2 seems to widen the gaps between the curves t = constant and to greatly 
increase the period when f = 1 and | = 2. 

It will be noticed that in diagram XII the numbers on the curves t = constant indicate the 

value of |* ■ 

(5) The chief effect of a decrease in f seems to be a slight change in curvature of the curves 
f = constant. 
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APPENDIX TO REPORT NO. 80. 

ABILITY OF THE PARACHUTE AND HELICOPTER. 
Note I. 

The expressions for H and K in terms of x 1} x it x„, y u y 2 , y 3 may be found as follows: 
Writing 

ft (v) = (v + Vi) (v + Vi) (y + y») - + y i) - W da + v J - ^(y + yO > 
/ 2 (y ) - acj (2/ + y 2 ) (y + Vi) +x 2 (y+ y s ) (y + yi) + x, (y+ yd (y + y 2 ) - aw„ 

we find that 
where 

(:x 1 + x 2 + x s ) 2 f 3 (y) =Ax 2 x 3 (y+y l ) +Bx i z l (y + y t ) + Cx^iy + y 3 ) . 

The coefficient of y on the right-hand side is the quantity denoted by //. 
Again writing 

fM-<M + W&)-K 

where a, I, and K are constants to be determined, we may find the value of K by substituting 
in the last equation a value of y which makes f s (y) =0. We thus find that 

Hiy + yJ^x^Bn-Cm), H(y+y 2 )=x 2 (Cl-An), H{y + y s )~x s {Am,-BT) 

= xJ J =z 2 <) = x 3 R 

H*K= xp&lE? -QB-RP-PQl 
Putting =0 in the identical relation connecting/^), f 2 (y), and/ a (y) we also find that 
I-PK(x 1 2 P + x 2 2 Q + x^B)^[H 2 (P + Q+B)-PQR]x l x 1 x 3 . 

HcilCG 

EPKUI+ixSP + ixfQ + ixfB) = x 1 x 2 x 3 (H+iP)(S+iQ)(H+iR) 
IPK^H 1 + (x*P + zfQ + x t >R)*J-z ] &c t V3* + P^iH 2 + Q*) \H*+R 2 )K 

In this relation the positive signs must be taken for the square roots so that K is positive 
when x v x 2! and x s are positive. 

Note II. 

In the case of a helicopter rising or falling vertically it may be sufficient to take into ac- 
count the fin action and gyroscopic action of the lifting screw. 
Writing 

<2 = o, E= X u , F= -iX r , J=iN„, K= N r +iN p 
the period equation becomes 

X(X+u)(X+i/-f is) +xQ\.w+i) = 0 

where 

g 2 N u N r N v 
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Routh's conditions indicate that for stability the quantities 
H=uyz?+ (y+u) 1 (wx+«y) — x(y + u), 

and 

K= ^TP - RuzH- (y +«)xy 2 2 3 ] 

must be positive. » 
On the other hand we find that 

so that the conditions E>o and T>o are equivalent. 

The effect of gyroscopic action on stability has been estimated for the case of an airplane 
in rectilinear flight and found to be small. 11 The value adopted for N v was N P — IQ where I 
is the moment of inertia of the propeller about the axis of y and 0 its angular velocity about 
this axis. 

With 7=150 pounds-/? and 0 = 2xX^|jj^= 125.8 radians per second, this gave a value of 

N p of about 15m for an airplane of mass m = 1,300 pounds. In the case of the helicopter, 0 is 
smaller than for the airplane propeller but I is very much greater if the diameter of the lifting 
screw is large. It seems likely, then, that the gyroscopic effect on stability will be greater 
than in the case of the airplane. 



'iL. Bairstow,B. SieiriE Jones, and A. W. H. Thompson, British Advisory Committea's Report. 1912-13, p. 166. ** 



